> 



A numerical study on the dimension of an 
extremely inhomogeneous matter distribution 



Institute) de Fisica, Universidade de Sao Paulo 
C.P.66.318, CEP 05315-970, Sao Paulo, Brazil 



in . 

O . Cecilia B. M. H. Chirenti 

<N 

< 

f"^ ■ Abstract 

We have developed an algorithm that numericaly computes the di- 
mension of an extremely inhomogeneous matter distribution, given by a 
discrete hierarchical metric. With our results it is possible to analise how 
the dimension of the matter density tends to d = 3 , as we consider larger 
samples. 

In a previous work HQ, we have presented a study on the hierarchical metric 
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dS 2 = dt 2 + g n dx 2 + g 22 dy 2 + g 33 dz 2 , (1) 

where the metric is defined on all integers, depending on their decomposition in 
' terms of powers of 2 as 
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922(y) = a(t) 2e , with y = 2 e+1 n + 2 e -l (2) 
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We obtained the following expression for the matter density, 

"S 1 la 2 

Too = P = 3-^-3 (k£ + £m + mk) 
87rG a z 

= p (t) (M + lm + mk) , (3) 

by means of the Einstein equations. Computing the Christoffel symbols and 
subsequently the curvature tensor for this metric requires some care, since we 
are not dealing with derivatives of functions, but differences of functions defined 
on a discrete space. 

We speculated whether such a matter distribution could be described by a 
fractal. As it turned out, from our preliminary analisis (see pQ), the dimension 
of the matter density tends slowly to d — 3. 
Considering the relation 

li m ^§1 = K , with N(r) = V — , (4) 

0<x,y,z<r ru 
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where K is some constant, and the constant d is the fractal dimension of the 
matter density 2; , it is easy to see that it implies 

\nN(r)=\nK + d\nr, (5) 

for large r. We have numerically computed N(r) and the plot In N(r) x lnr, in 
figure QJ, showing that equation (J5J, although valid at very large r), is a good 
approximation for the behavior of the data. 
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Figure 1: In this graphic we have In N(r)x lnr. The equation of the straight line f(x) = ax+b 
has a = 3.00374 ± 0.000025 , b = 1.06288 ± 0.000225 . 

Numerical fits for these points, taking increasingly larger samples, have given 
results for d that approach d = 3, as shown in figure (0- 
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Figure 2: This graphic shows (d — 3) X lnr, obtained for r values up to 20.000. The 
curve is a parabola fitted to the points by Gnuplot, with equation h(x) = ax 2 + bx + c: a = 
0.00102 ± 0.00014 , b = -0.02134 ± 0.00220 , c = 0.11576 ± 0.008441 . 

As a conclusion, we can say, now, that the dimension of the matter density 
that generates metric (QJ tends to d = 3 as shown in figure (J5J) . 
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